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Abstract. Boundary conditions in quasiclassical theory of superconductivity are of 
crucial importance for describing proximity effects in heterostructures between different 
materials. Although they have been derived for the ballistic case in full generality, 
corresponding boundary conditions for the diffusive limit, described by Usadel theory, 
have been lacking for interfaces involving strongly spin-polarized materials, such as e.g. 
half-metallic ferromagnets. Given the current intense research in the emerging field of 
superconducting spintronics, the formulation of appropriate boundary conditions for 
the Usadel theory of diffusive superconductors in contact with strongly spin-polarized 
ferromagnets for arbitrary transmission probability and arbitrary spin-dependent 
interface scattering phases has been a burning open question. Here we close this 
gap and derive the full boundary conditions for quasiclassical Green functions in the 
diffusive limit, valid for any value of spin polarization, transmission probability, and 
spin-mixing angles (spin-dependent scattering phase shifts). It allows also for complex 
spin textures across the interface and for channel off-diagonal scattering (a necessary 
ingredient when the numbers of channels on the two sides of the interface differ). As 
an example we derive expressions for the proximity effect in diffusive systems involving 
half-metallic ferromagnets. In a superconductor/half-metal/superconductor Josephson 
junction we find ((ig-junction behavior under certain interface conditions. 
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1. Introduction 

Hybrid structures containing superconducting (S) and ferromagnetic (F) materials 
became a focus of nanoelectronic research because of their relevance for spintronics 
applications as well as their potential impact on fundamental research mw- Examples 
of successful developments include the discoveries of the vr-junction H E] in S/F/S 
Josephson devices iiEi, of odd-frequency superconductivity [8] in S/F heterostructures 
PHO], and of the indirect Josephson effect in S/half-metal/S junctions [IT[[T2]. Other 
recent topics of interest include the study of Majorana fermions at interfaces between 
superconductors and topological insulators ra and at edges in superfluid ^He [HI |T2] , 
and the appearance of pure spin supercurrents in topological superconductors [I6] , and 
in S/FI-F-FI devices as a result of geometric phases [T7j . 

The central subject in many of these studies is to understand how in the case 
of a superconductor coupled to a ferromagnetic material superconducting correlations 
penetrate into the ferromagnet, and how magnetic correlations penetrate into the 
superconductor [EliiniEniEllEaEa]. A powerful method to treat such problems is the 
quasiclassical theory of superconductivity developed by Larkin and Ovchinnikov and by 
Eilenberger |2ll [25] . Within this theory [26l EH EHl ES ED] the quasiparticle motion is 
treated on a classical level, whereas the particle-hole and the spin degrees of freedom 
are treated quantum mechanically. The transport equation, which is a hrst order matrix 
differential equation for the quasiclassical propagator, must be supplemented by physical 
boundary conditions in order to obtain a unique solution. 

Whereas for the full microscopic Green functions, the Gor’kov Green functions 
[DT] , such boundary conditions can be readily formulated {e.g. in terms of interface 
scattering matrices or in terms of transfer matrices), this is a considerably more difficult 
task for quasiclassical Green functions. In quasiclassical theory only the information 
about the envelope functions of Bloch waves is retained, information about the phases 
of the waves is however missing. Such envelope amplitudes can show jumps at interfaces, 
and one complex task is to calculate these jumps without knowing the full microscopic 
Green functions near the interface. Gorrespondingly, there is a long history of deriving 
boundary conditions for quasiclassical propagators, both for the Eilenberger equations, 
and their diffusive limit, the Usadel equations |32j . 

For ballistic transport, described by the Eilenberger equations, such boundary 
conditions were hrst formulated for spin-inactive interfaces in pioneering work by 
Shelankov and by Zaitsev m l35] . who showed the non-trivial fact that these jumps 
can be calculated using only the envelope functions. More general formulations were 
proposed subsequently [36l [371 EHl EH], including a formulation in terms of interface 
scattering matrices by Millis, Rainer, and Sauls [39] • All these formulations were 
implicit in terms of non-linear matrix equations, and problems arose in numerical 
implementations due to spurious (unphysical) additional solutions which must be 
eliminated. Progress was made with the help of Shelankov’s projector formalism 
ra. allowing for explicit formulations of boundary conditions in both equilibrium 
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[m 02], 03] and non-equilibrium |12] situations. Further generalizations included spin- 
active interfaces, formulated for equilibrium [H] and for non-equilibrium [l5], and 
interfaces with diffusive scattering characteristics |16|. An alternative formulation in 
terms of quantum mechanical t-matrices [17| proved also fruitful [HI 09] [HI jJO] (50] [51] . 
The latest formulation, in terms of interface scattering matrices, is able to include non- 
equilibrium phenomena, interfaces and materials with weak or strong spin polarization, 
multi-band systems, as well as disordered systems |52j . 

For the diffusive limit a set of second order matrix differential equations has been 
derived by Usadel [32] • In contrast to the ballistic case, where boundary conditions have 
been formulated for a wide set of applications, boundary conditions for the diffusive limit 
have been formulated so far only in certain limiting cases. The first formulation is by 
Kupriyanov and Lukichev, appropriate for the tunneling limit ]H] • This was generalized 
to arbitrary transmission by Nazarov A major advance was done by Cottet et al 
in formulating boundary conditions for Usadel equations appropriate for spin-polarized 
interfaces [55]. These boundary conditions are valid in the limit of small transmission, 
spin polarization, and spin-dependent scattering phase shifts (this term is often used 
interchangeably with “spin-mixing angles” ]SH])- Subsequent formulations allowed for 
arbitrary spin polarization, although being restricted to small transmission and spin- 
dependent scattering [SH EH] ES]. In Ref. [H] the authors present “heuristically” 
deduced boundary conditions, which coincide with the ones used in Refs. [57] EH]. 

Here we not only present the full derivation of the specific boundary conditions 
used in Refs. [571EHIE9], but go further and give a full solution of the problem. With 
this, the long-standing problem of how to generalize Nazarov’s formula for arbitrary 
transmission probability [51] to the case of spin-polarized systems with arbitrary spin 
polarization and arbitrary spin dependent scattering phases is solved. Our boundary 
conditions are general enough to allow for non-equilibrium situations within Keldysh 
formalism, as well as for complex interface spin textures. We reproduce as limiting 
cases all previously known formulations. 

2. Transport Equations 

The central quantity in quasiclassical theory of superconductivity [211 [ 21 ] is the 
quasiclassical Green function (“propagator”) ^(p^,R,E,f). It describes quasiparticles 
with energy E (measured from the Fermi level) and momentum p^? moving along 
classical trajectories with direction given by the Fermi velocity v^(pj?) in external 
potentials and self-consistent fields that are modulated by the slow spatial (R) and 
time {t) coordinates [2HI [27] [2H] . The quasiclassical Green function is a functional of self¬ 
energies S(pi?, R, E, f), which in general include molecular fields, the superconducting 
order parameter A(pi?,R, f), impurity scattering, and the external potentials. The 
quantum mechanical degrees of freedom of the quasiparticles show up in the matrix 
structure of the quasiclassical propagator and the self-energies. It is convenient to 
formulate the theory using 2x2 matrices in Keldysh space [HD] (denoted by a “check” 
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accent), the elements of which in turn are 2x2 Nambu-Gor’kov matrices [6T] in 
particle-hole (denoted by a “hat” accent) space. The structure of the propagators and 
self-energies in Keldysh-space is 


9 = 


9^ 9 
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0 


s = 


9 


kel 


0 



(la) 


where the superscripts R, A, and K refer to retarded, advanced, and Keldysh 
components, respectively, and with the particle-hole space structure 
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for self-energies. For spin-degenerate trajectories {i.e. in systems with weak or no 
spin-polarization) the elements of the 2x2 Nambu-Gor’kov matrices are 2x2 matrices 
in spin space, e.g. with a,b E {t, i}, and similarly for others. In strongly 

spin-polarized ferromagnets the elements of the 2x2 Nambu-Gor’kov matrices are spin- 
scalar (due to very fast spin-dephasing in a strong exchange field), and the system must 
be described within the preferred quantization direction given by the internal exchange 
field. The terms “weak” and “strong” refer to the spin-splitting of the energy bands 
being comparable to the superconducting gap or to the band width, respectively. In 
writing Eqs. (la)-(lc) we used general symmetries, which are accounted for by the 
“tilde” 


operation, 

X(pF,R,F;,f) = X(-pj7’,R, 


-E,ty 


( 2 ) 


Retarded (advanced) functions can be analytically continued into the upper (lower) 
complex energy half plane, in which case the relation is modified to X(pi7’,R, E,t) = 
X{—pF,'R,,—E*,t)* with complex E. 

The quasiclassical Green functions satisfy the Eilenberger-Larkin-Ovchinnikov 
transport equation and normalization condition 

[Ef 3 -t,g]^ + ih-VF-Vg = 6, g o g = (3) 

The non-commutative product o combines matrix multiplication with a convolution over 
the internal energy-time variables in Wigner coordinate representation, 

{AoB){E,t) =eA9^^?-9t^E)A{ER)B{ER), (4) 

and fs = fsl, where % is a Pauli matrix in particle-hole space. Here and below, 
[A, B]^ = A o B — B o A. The operation V acts on the variable R. 


^ For the definitions of all Green functions in this paper we use a basis of fermion field operators in 
Nambu ® spin-space as 'I'(r,t) = t),'0^(r, t),t)b V’t(''R)^]^ ■ 
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The functional dependence of the quasiclassical propagator on the self-energies is 
given in the form of self-consistency conditions. For instance, for a weak-coupling, s- 
wave order parameter the condition reads 

/ Ec JJ 7 

^^{NF{pF)fJ{pF,R,E,t))p,, (5) 

'Ec 

where Vg is the s-wave part of the singlet pairing interaction, Np is the density of states 
per spin at the Fermi level, is spin-singlet part of the the Keldysh component f^, 
and ( )p^ denotes averaging over the Fermi surface. The cut-off energy Ec is to be 
eliminated in favor of the superconducting transition temperature in the usual manner. 

When the quasiclassical Green function has been determined, physical quantities 
of interest can be calculated. For example, the current density at position R and time 
t reads (with e < 0 the electron charge) 

/ OO 

—TT{Np{pF)vp{pF)fsg^{pF, R, E, f))p^. (6) 

-OO Svri 

The symbol Tr denotes a trace over the 2x2 particle-hole space as well as over 2x2 spin 
space in the case of spin-degenerate trajectories. 

In the dirty (diffusive) limit, strong scattering by non-magnetic impurities 
effectively averages the quasiclassical propagator over momentum directions. The Green 
function may then be expanded in the small parameter k-QTcT/h (r is the momentum 
relaxation time) following the standard procedure [32l |33] 

g{pF, R, E, t) ^ (^(R, E, t) + g^^\pF, R, E, t) (7) 

where the magnitude of g^^'^ is small compared to that of G. The impurity self-energy is 
related to an (in general anisotropic) lifetime function r(p'^,pi7’) |33]. Substituting 
([^ into ([^, multiplying with iV^(p)7’)vi7’j(p);.)r(p);., PiT’), averaging over momentum 
directions, considering that S'r/h is small, where S' is the self-energy reduced by the 
contribution due to non-magnetic impurity scattering, and using G o G = —vr^l and 
G o g^L _|- gW o G = 6, one obtains (we suppress here the arguments R, E, t) 

[NF{pF)^F,j{VF)g^^\vF)) = Np ^ -^G o VfcG, (8) 

k 

where Np = {Nf{pf))pp is the local density of states per spin at the Fermi level, 
Vfc = d/dRk, the summation is over k G {x,?/, x}, and 

Djk = ^((NF{p'F)^F,j{p'F)^{p'F^PF)yF,k{PF)NF{pF)) ( 9 ) 

JVp, \ \ ' PF' Vf 

is the diffusion constant tensor. For isotropic systems, Djk = D5jk- The Usadel Green 
function G obeys the following transport equation and normalization condition [32], 

[Ff3i-So,G], + 5^^V,(GoVfcG)=0, GoG = -f% (10) 

jk 

where So = {Nf{pf)^'{pf))pf/E^f- The Usadel propagator G is a functional of Sq. 
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The structures of G and Sq are the same as in Eqs. (la)-(lc) (with G replacing g 
and So replacing S). Eq. ([^ is replaced by 

X(R, E, t) = X(R, -E, ty. (11) 


The current density for diffusive systems is obtained from Eqs. (|^ and (|^, and is given 
by 


ji(R,t) = -e^ f 
1^ ^ — 


—TTNpD,kfyG{R,E,t)oVkG{R,E,t)f. ( 12 ) 


A vector potential A(R, t) enters in a gauge invariant manner by replacing the spatial 
derivative operators in all expressions by (see e.g. [331IS2]) 


^ a o X = -1 


h 


(13) 


Finally, the case of a strongly spin-polarized itinerant ferromagnet with 
superconducting correlations {e.g. due to the proximity effect when in contact with 
a superconductor) can be treated by quasiclassical theory as well [HI [201 EO]- In this 
case, when the spin-splitting of the energy bands is comparable to the band width of 
the two spin bands, there exist two well separated fully spin-polarized Fermi surfaces in 
the system, and the length scale associated with h/|pF'|- — is much shorter than the 
coherence length scale in the ferromagnet. Equal-spin correlations stay still coherent 
over long distance in such a system; and correlations are, however, incoherent 
and thus negligible within quasiclassical approximation. Fermi velocity, density of 
states, diffusion constant tensor, and coherence length all become spin-dependent. The 
quasiclassical propagator is then spin-scalar for each trajectory, with either all elements 
or all elements j,j, depending on the spin Fermi surface the trajectory corresponds 
to. Eilenberger equation and Usadel equation have the same form as before for each 
separate spin band. The spin-resolved current densities are given in the ballistic case by 

/ OO 

—Tr(XFtVFtr3^n)pFt’ 


and in the diffusive case by 


dE 


Jfct = 


k ' 




K 


(15) 


and analogously for spin down. 

For heterostructures, the above equations must be supplemented with boundary 
conditions at the interfaces. A practical formulation of boundary conditions for diffusive 
systems valid for arbitrary transmission and spin polarization is the goal of this paper. 


3. Boundary Conditions 

3.1. Interface Scattering Matrix 

We formulate boundary conditions at an interface in terms of the normal-state interface 
scattering matrix S [3Sl EH ES], connecting incoming with outgoing Bloch waves on 
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either side of the interface with each other. We use the notation 



where 1 and 2 refer to the two sides of the interface, and the subscript label indicates 
that the 2x2 matrix structure refers to reflection and transmission amplitudes at an 
interface. The components Sij are matrices in particle-hole space as well as in scattering 
channel space {e.g. scattering channels for ballistic transport would be parameterized 
by the Fermi momenta of incoming and outgoing Bloch waves). Each element in 2x2 
particle hole space is in turn a matrix in combined spin and channel space, i.e. the 
number of incoming directions (assumed to be equal to the number of outgoing directions 
due to particle conservation) gives the dimension in channel space. The dimension in 
spin space is for spin-degenerate channels 2 and for spin-scalar channels 1. 

If time-reversal symmetry is preserved, Kramers degeneracy requires that each 
element of the scattering matrix has a 2x2 spin (or more general: pseudo-spin) structure 
(as it connects doubly degenerate scattering channels on either side of the interface). 
For spin-polarized interfaces {e.g. ferromagnetic or with Rashba spin-orbit coupling) the 
scattering matrix is not spin-degenerate. However if the splitting of the spin-degeneracy 
is on the energy scale of the superconducting gap, it can be neglected within the precision 
of quasiclassical theory of superconductivity. On the other hand, if the lifting of the 
spin-degeneracy of energy bands is comparable to the Fermi energy, the degeneracy of 
the scattering channels must be lifted as well in order to achieve consistency within 
quasiclassical theory. For dehniteness, we denote the dependence on the scattering 
channels by indices n, n'\ 


[Sa/3]nn', (17) 

even for the ballistic case for which [Sa/^jnn' = Sa^(pi?^„, ki?_„/). 

As shown in Appendix A and [Appendix B the scattering matrix for an interface 
can be written in polar decomposition in full generality as 


Vi - ccf c A / 5 0 A 
-Vl - cw j I 0 5 j 

/x ^ Z X 

with unitary matrices S and <S, and a transmission matrix C. All are matrices in 
particle-hole space, scattering channel space, and possibly (pseudo-)spin space. The 
above decomposition divides the scattering matrix into a Hermitian part and a unitary 
part. From this decomposition, we can dehne the auxiliary scattering matrix 





5 0 
0 5 


(19) 


which retains all the phase information during reflection on both sides of the interface, 
and has zero transmission components. The decomposition is uniquely dehned when 
there are no zero-reflection singular values (we will assume here that always a small 
non-zero reflection takes place for each transmission channel; perfectly transmitting 
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channels can always be treated separately as the corresponding bonndary conditions 
are trivial). For the matrix C we introdnce the parameterization 


C = (l + tt^) ^ 2t, 


( 20 ) 


see 


Appendix C) which is nniqnely defined when all singnlar valnes of t are in the 


interval [0,1] (which is required in order to ensure non-negative reflection singular 
values). We define for notational simplification “hopping amplitude” matrices 

7rri2 = ts, 71T21 = t^S, (21) 

as well as unitary matrices 

5i = 5, ^2 = 5. (22) 

In terms of those, obviously the relation 

Taa = (23) 

holds, where {a, a) G {(1,2), (2,1)}, and the labels 1 and 2 refer to the respective 
sides of the interface. Here, and below, the Hermitian conjugate operation involves a 
transposition in channel indices. The particle-hole structures of the surface scattering 
matrix and the hopping amplitude are given by. 





'T~rvn/ — 


'^aa 
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(24) 


ph 


with 


[‘S'ojnn' [haajnn' [haalfifi '5 (2^) 

where n and n' denote mutually conjugated channels, e.g. defined by pF,n' = —k^.n' 
and kp^n = ~PF,n- Finally, the Keldysh structure of these quantities is 
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^rvr\ 



0 \ 


[ 0 

) 

kel \ 

f f^- 
1 ' 0 . 0 . 

0 ^ 


1 •> 


le, ^ ^ 



(26) 

(27) 


(the additional Hermitian conjugate in these equations is due to the fact that advanced 
Green functions have the roles of “incoming” and “outgoing” momentum directions 
interchanged compared to retarded Green functions; this is similar to the additional 
Hermitian conjugate appearing for hole components in particle-hole space). Thus, the 
Keldysh matrix structure for Sa and faa is trivial (proportional to unit matrix). The 
full normal-state scattering matrix is diagonal in particle-hole and in Keldysh space, 
with reflection components 

S«a = (1 + (1 - Sa, (28) 

and with transmission components 

Seta (1 “h TT 27rTaa- 


(29) 
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Figure 1. (a): Illustration of notation used in this paper, (b) and (c): Structure of 

boundary condition with transfer matrices M in (b), and with scattering matrices S 
in (c) (yellow). “Drone” amplitudes in the propagators (orange fields) connect in (b) 
incoming (i) and outgoing (o) momentum directions, and in (c) the two sides, a and 
a, of the interface. To obtain quasiclassical boundary conditions. Drone amplitudes in 
(b) and (c) must be eliminated. In this paper we use formulation (c). To connect to 
the notation in the main text, = g\ g'‘f^ = g\ g°ff = g°, and g°ff = g°. 


Note that Taa connects incoming with outgoing Bloch waves per dehnition (as the 
scattering matrix does). 

We will formulate the theory such that all equations are valid on either side of the 
interface. This allows us to drop the indices a, a for simplicity of notation by randomly 
choosing one side of the interface, and denoting quantities on the other side of the 
interface by underline. In particular, we will use 


Set — 5*, Sa — S_, Taa — L ; Tcfa — T. 

9a = 9i 9a = fr Ga = G, Ga = G, 
and so forth [see figure [^a)]. Also, from Eq. (23) we have f = Sf^^. 


(30) 


3.2. General Boundary Conditions for diffusive systems 

One main problem with boundary conditions for quasiclassical propagators is illustrated 
in hgure (b) and (c). In previous treatments [39l |Ml |55] the starting point was a 
transfer matrix description, see hgure (b), which however required the elimination of 
so-called “Drone amplitudes”, which are propagators that mix incoming with outgoing 
directions. Here, we will employ a scattering matrix description, see hgure [^(c), which. 
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on the other hand, reqnires a similar elimination of Drone amplitndes, this time being 
propagators mixing the two sides of the interface. However, for an impenetrable interface 
this latter problem does not arise, a fact we will exploit. 

The strategy to derive the needed bonndary conditions is to apply a three-step 
procednre. In the first step, the problem of an impenetrable interface with the anxiliary 
scattering matrix defined in Eq. (19) is solved on each side of the interface m- For 


this step, the ballistic solntions for the envelope fnnctions for the Gor’kov propagators 
close to the interfaces shonld be expressed by the solntions G of the Usadel eqnation. 
In a second step, these ballistic solntions (anxiliary propagators) are used in order to 
find the full ballistic solutions for finite transmission by utilizing a f-matrix technique 
[IHllIIlEollso]. In the third, and final, step the matrix current will be derived from the 
ballistic solutions, which then enters the boundary conditions for the Usadel equations. 
We will present explicit solutions for all three steps, such that the procedure describes 
effectively boundary conditions for the solutions of Usadel equations on either side of 
the interface. 

We use for the auxiliary propagators the notation ^q, ^q, g° and , where the upper 
index denotes the direction of the Fermi velocity. Incoming momenta (index i) are those 
with a Fermi velocity pointing towards the interface, and outgoing momenta (index o) 
are those with a Fermi velocity pointing away from the interface. 


3.2.1. Solution for impenetrable interface: We solve first for the auxiliary ballistic 
propagators fulhlling the impenetrable boundary conditions 

d 


r« = SgiSl, g = SgS', (31) 

implying matrix multiplication in the combined [Keldysh] x [particle-hole] x [combined 
scattering-channel and spin] space. For diffusive banks, it is necessary to connect the 
ballistic propagators gQ° with the isotropic solutions of the Usadel equation, G. The 
ballistic propagators and g‘^^°, which characterize electronic correlations next to the 
scattering barrier, depend on the electronic momentum. However, in the diffusive case, 
impurity scattering leads to momentum isotropization away from the scattering barrier. 
This process occurs in isotropization zones with a thickness corresponding to a few 
times the inelastic mean free path of the materials, see figure (a). This scale is itself 
much smaller than the scale on which the isotropic diffusive Green functions evolve 
in the bulk of the materials, in the framework of the Usadel equations. Indeed, the 
Usadel equations involve a superconducting coherence length, which is typically much 
larger than the elastic mean free path. Therefore, in order to describe disordered hybrid 
structures with Usadel equations, suitable boundary conditions should be expressed in 
terms of the values of the isotropic Green functions G and G right at the beginning of the 


isotropization zones. To obtain such boundary conditions from Eq. (31), it is necessary 


to express the propagators ^q’ and in terms of G and G. This can be done by 
studying the spatial dependence of the Gor’kov Green functions (or full Green functions 
without the quasiclassical approximation) in the isotropization zones (see Refs. [HI EH] 
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for details). Using the fact that the dynamics of electrons is dominated by impurity 

scattering in these zones, one can express the Gor’kov Green functions in terms of , 

G and G. Then, an elimination of unphysical solutions imposes the conditions 0 

{G — ittI) o + ittI) = 0, {gh — ivri) o (G + ittI) = 0 (32a) 

{G + ivri) o — ivri) = 0, (^g + ° {G — i^ri) = 0 (326) 

and similarly for G and g^^°. From this one obtains the identity ^ = —vr^l for 

the anticommutator This allows to solve after some straightforward algebra for 

using Eq. (31), and using the abbreviations 


G' = ^ (S'GS 


G), G" = (SOS' - G), (33) 

(both are matrices depending via S on the scattering channel index) leading to [55] 

— ittI = (1 — G o G')“^ o (G — ittI), (34a) 

^« + i7ri = (l-GoG")“^o(G + i7ri) (346) 

(here and below the inverse is defined with respect to the o-product), which, using 


identities like G' o G' = ^ 

can be written also as 

gl + ivri = (G + ittI) o (1 
(Jq — ivrl = (G — ittI) o (1 


(with {A, B}^ = A o B + B o A), alternatively 


G'o 

G"c 


G)-\ 

<G)-\ 


(34c) 
(34 (i) 


Similarly equations hold for G and g^^° in terms of the scattering matrix S_. Introducing 
these solutions into Eqs. (32a)-(326) shows readily that the latter are fulhlled. We note 
that the relation ^q° o ^g’" = —vr^l follows from G o G = —tt^I and SS'^ = S'^S = 1. It is 
also important to notice that whereas G is proportional to the unit matrix in channel 
space due to their isotropic nature [5^, S, and consequently G', G", and ^g°, are in 
general non-trivial matrices in channel space. Eqs. (34a)-(346), or alternatively (34c)- 

in terms of the diffusive Green 
re difference ^g — ^g in a more explicit manner, using the 


(34d), together with Eq. (33) determine uniquely g^ 
function G 


We can rewrite t 

abbreviations S' = G o G' and S" = G" o G, leading to 


9o 


= (i — (5') ^ o [(G — ivri) o 5" — 5' o (G — ivri)] o (1 — 5") (35) 


3.2.2. Solution for finite transmission: The second step follows Refs. HHEQ]. Once 
the auxiliary propagators are obtained, the full propagators can be obtained directly, 
without further solving the transport equation, in the following way. We solve t-matrix 
equations resulting from the transmission parameters f, for incoming and outgoing 
directions, which according to a procedure analogous to the one discussed in Ref. mm 
take the form, 

t" = i o (l + ^0 ° t") ) ^ io ° + ^0 ° i°) . (36) 

Using the symmetry Eq. ( [2^ , the t-matrices for incoming and outgoing directions can 
be related through 

= Si' Si 


(37) 
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Using the short notation 


4 


9l = rgf^T\ 9\=t^^D (38) 

we solve formally Eqs. ( |3^ for f’°: 

( 39 ) 

The full propagators, fulhlling the desired boundary conditions at the interface, can 
now be easily calculated. For incoming and outgoing directions they are obtained from 

[nlEo] 


. -1 


9 " = 9o+ {9o + iTTl) o ? o (^0 - ) (40a) 

g'' = 9o + (go - iTri) o t° o (g^ + iTri) . (406) 

Noticing that (^q° + ivri) o (g'^^° — ivri) = 0, and (^q° — ivri) o (g'^^° + ivri) = 0, as well as 
identities like ^q°o (^Q° + i7rl) = i7rlo(^Q° + i7rl) etc, it is obvious that the normalization 
gho o gi,o _ holds. Using the same identities, we obtain the alternative to Eqs. 


(40a)-(406) expressions 

g' =9o + {9o + iTTl) o [t\ gf\o = 90 - [t\ gflo o {g}^ - ivri), 
9° = 9o + {9o - iTTl) o [i°, go]o =90- [t°, go]o o (^o + i^rl). 


(40c) 
(40 d) 


Equations (40a)-(406), or alternatively, ( |40c )-(40d), in conjunction with Eqs. (38)-(39), 
solve the problem of Ending the ballistic solutions for finite transmission. We are now 
ready for the last step, to relate these solutions to the matrix current which enters in 
the expression for boundary conditions for G and G. 


3.2.3. Matrix current and boundary conditions for diffusive propagators: We now turn 
to the third, final, step. As shown in Refs. [541 155] , the boundary conditions for 
quasiclassical isotropic Green functions can be obtained from the conservation of the 
matrix current X in the isotropization zones surrounding the scattering barrier. This 
quantity contains physical information on the flows of charge, spin and electron-hole 
coherence in a structure. We refer the reader to Refs. [HU [55] for the general definition 
of X in terms of the Gor’kov Green functions. Using this definition, one can verify that X 
is spatially conserved along the entire isotropization zones. Then, one can express X next 
to the scattering barrier in terms of the propagators g^'° and g'''°, and at the beginning 
of the isotropization zones in terms of G and G, see Fig. (a). The conservation of 
the matrix current provides an equality between the two expressions. Since g^'° can be 
expressed in terms of go° and g‘^^°, and these in terms of the G and G, this gives the 
desired boundary conditions. Following Ref. [50], after some straightforward algebra 
we obtain 


9o]o = (1 - o g°) ^ [gl, gf\^ (1 - ° 9?) 




Using relations (31) and (37) above, we find 


(41) 


9 


[^0 + (^0 + ittI) o i° o (g° - ivri)] S, 


(42) 
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which allows to derive the following relation 

X'^gO-Sg^S^= -27:i[i°G9l]o. (43) 

For calculating the charge current density in a given structure, it is sufficient to know 
X', because the matrices S and drop out of the trace as they commute with the fa 
matrix in particle-hole space. 

Finally we relate the obtained propagators to the matrix current X, 

X = g°-g^=X'+X" (44) 

with 


i" = Sg‘S' - g‘. 


(45) 


We remind the reader here that X has a matrix structure in Keldysh space, in 
particle-hole space, and in combined scattering-channel and spin space. In terms of X the 
boundary condition results then from Eq. (|^ and from the matrix current conservation 
in the isotropization regions ig 



71=1 


(46) 


where 2 ; is the coordinate along the interface normal {away from the interface), n 
is a scattering channel index {M channels, spin-degenerate channels count as one), 
cr = e^NpD refers to the conductivity per spin, A is the surface area of the contact, 
and Qq is the quantum of conductance, Qq = jh. The number of scattering channels 
is expressed in terms of the projection of the Fermi surfaces on the contact plane, 
by W = Ap^zAj{2 'kY. For isotropic Fermi surfaces Ap^z = In general, 

1 r d‘^k\\ 

where Kk\\ is the momentum component parallel to the interface. 



4. Special Cases 


4-1. Spin-scalar and channel-diagonal case 

The transition to the diffusive Green functions is trivial for the case of S' = 1, as then 


9o — 9q — If start from Eq. (41) in conjunction with (38), we obtain in the case 
of a spin-scalar and channel-diagonal matrix f„„ with the notation G = —ittG 




ITT 


E 


4r„[G,G] 


4 + r„({G,G}-2) 6, 


^^-^G o Tg 


dz 


with a = e^NpD and 


r„ = 


47r^ It 

^* \ ‘ n. 


(1 7r2|r„„|2) 

This reproduces Nazarov’s boundary condition 


(48 a) 


( 486 ) 
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For the case of zero transmission, f = 0, we can find a closed solntion if we assume 
that we can hnd a spin-diagonal basis for all reflection channels. For a channel-diagonal 
scattering matrix we write Snn = with k = diag {mcf, mcf*}, where = 1 

(leading to = 1). In this case we have g^’° = gl’^. We use Eq. (35), which 
straightforwardly leads to 

2 Yin 


ITT 




- isin?9„ - - 

1 — ---(G/xCt — K) -|- 


sin^ ^ 


{GkGk-l) 


X 


19 

c r, 


-isin'd„[/t, G] -|- sin^ -^[kGk, G] 


X 


1 sin 'd. 


gj]^2 iLi. 

-{GkG -k) + -^^{kGkG - i) 


-1 


(49) 


(where we remind that G^ = 1). Note that (pn drops out, only the spin mixing angle 


dn matters. Eq. (49) generalizes the results of Ref. |55] to arbitrary spin-dependent 
reflection phases. Further below we will give a physical interpretation of the leading 
order terms arising in an expansion for small dn- 


4-3. Exact series expansions 

We now provide explicit series expansions for all quantities which will be useful for 
deriving formulas for various limiting cases. We start with writing the scattering matrix 
as ^ with hermitian K due to unitarity of S, i.e. K = KY Then we use an 

expansion formula for Lie brackets in order to obtain the series expansion 


S^GS = e 




>-i 


m=0 


ml 


[K^G] 


(50) 


with the dehnitions [iF T [K ™', ^ G]] and [iF (• G] = G. With this we obtain 


from Eq. (33) 


m=l 


G" = 




2'7r^ —' m\ 

m=l 


( 61 ) 


which are very useful if K has a small pre-factor. Note also the identity Go |^iF, G] oG = 
TT^ [iF, G] . Furthermore, from Eqs. (^34cH34d) we hnd 


% = 


G -|- (G -|- ittI) o E(G' ° G)' 


1=1 

oo 


^« = G + (G-ivri)o^(G"oG)'. 


1=1 


From Eq. (41), and using Eqs. ( |31| ), ( |37[ ), we derive 

OO 

^ o ] o = E ° ° ^ o ]° ° *^^0 ° 


(52 a) 
(526) 

(53 a) 


/c,n=0 
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= J2^9\°9o)’'°[9l9o],°{9o°9\T, (536) 

k,n=0 

which is useful if the transmission amplitudes f entering into gi° are small. Finally, we 


obtain from Eqs. (43) and (45) 

OO ™ 


m=l 


Here, is obtained from 

OO 

+ ivri = (G + ittI) o o G)' o (i + [t*,^o]o) 


(54) 


(55) 


1=0 


4-4■ Boundary condition for spin-polarized surface to third order in spin-mixing angles 

We hrst treat the case when f’° = 0, for example the case where one side of the junction 
is a ferromagnetic insulator (FI). Then 


TTl 771 


( 66 ) 


m=l 


m,/=l 


To third order we have X = and the derivation in Appendix D leads to 

(57a) 

(576) 


x(^) = 1 [k, G] , x(2) = -^ [iFGi^, G] 


X(3) = -^[iF3G]-^[iF,Go[iF?G]oG] 


24 ’ J 87r2 

For the special case of channel diagonal knn = with ik' = 1^ which follows also from 


directly expanding Eq. (49), we reproduce the results from Ref. [55] (G = —fvrG), 


2En^ 


( 1 ) 


ITT 


= -1 (En'^n) [k, G] , 


2 En En 


2 

n '' n 


ITT 


2E„i: 


(3) 

7177 


= —1- 


. /-^77 77 


- [k,G] - [kGkoGk^G]^ 


[kGk, G]^ (58a) 


(586) 


ITT 16 \3 

Note that the hrst order term ~ [k, G] accounts for the effective exchange held induced 
inside the superconductor by the spin-mixing, whereas the term ~ [kGk, G] produces 
a pair breaking ehect similar to that of paramagnetic impurities [66|. This second term 
occurs only at second order in {in because it requires multiple scattering at the S/FI 
interface, which together with random scattering in the dihusive superconductor leads 
to a magnetic disorder ehect. 

4-5. Boundary condition for spin-polarized interface to second order in spin-mixing 
angles and transmission probability 

We now allow for hnite transmission, and concentrate on the matrix current to second 
order in the quantities iF, and ki°- We need to take care of the scattering phases 
during transmission events. For this, we dehne 

f = ShoS\ f = S^roSX 


(59) 
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We note that Eq. (23), or f = Sf^S_, results into 


T-o = 


(60) 


Thus, the Tq and fg are the appropriate transmission amplitudes, with transmission 
spin-mixing phases removed. We further define 


Gi = tqGtq. (61) 

We expand f up to first order in K and K, 

f = fg ^ {kfo fok) + (62) 


and obtain X = X*^^) -I- X*^^) from a systematic expansion to second order in K, k and 
Gi, as shown in Appendix E, leading to one of the main results of this paper 


XW= -27ri[Gi,G]^ + i[i^,G], 

x(2) = _ 27ri [Gi O G O Gi, G] ^ ^ [KGk, G] ^ 

+1 [Gi o Gk + i^G o Gi + fgG o [k G] fg^ G 


(63 a) 


(636) 


These relations generalize the results of Ref. [55] for the case of arbitrary spin 
polarization, and are valid even when K, k and r have different spin quantization 
axes, i.e. cannot be diagonalized simultaneously. 

Using the notation G = —ittG and 27rfo = T, we can rewrite the result in leading 
order in the quantities K, k and the transmission probability ( fft) as 

2X(^) , , , . . 

— = [TGrt-2iA',G]^, (64a) 

and for the next to leading order 

2T(2) 1. 

= - - [TGkoGoTGk,G]^ + [KGK,G]^ 

+ ^[TGkoGk + kGofGk + TGo [X,G]f^G]^. (646) 

These equations are still fully general with respect to the magnetic (spin) structure, 
and allow for channel off-diagonal scattering as well as different numbers of channels on 
the two sides of the interface. Note that T, K, and k are matrices in channel space, 
whereas G and G are proportional to the unit matrix in channel space. Whereas K, 
and k are square matrices, T in general can be a rectangular matrix (when the number 
of channels on the two sides of the interface differ). 


4-6. Boundary conditions for channel-independent spin quantization direction 

As an application, we assume next that each of the quantities K, k and f can be 
spin-diagonalized simultaneously for all channels, with spin quantization directions m', 
fn, and rh for K, k, respectively. We also use that G and G are proportional 
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Tnn' 1 + ^^nn' rfl' ■ & = Knn', 1 + fu ■ & = KlV, 


a = (tI, a = 


to the unit matrix in channel space, as they are isotropic [55], and we assume that the 
number of channels on both sides of the interface are equal. We dehne 

To,nz 1 + Ti^ni rfl - a = tni, (65a) 

—- ^ 

^ ^ ^1 ^ ■. 

^ ^^1 , K = m ■ a, K = m ■ a, k = m ■ a (65c) 

/ ph 

with = (m')^ = {fn)‘^ = 1, i.e. k‘^ = {k'^ = (k')^ = 1, and introduce the 
transmission probability Tni and the spin polarization Vni as 

Tnl (l + Vnirha) = fnl[fnl\'. (66) 

We write for Tq, nz and allowing for some spin-scalar phases V’nz, 

(67) 


pjp2 ' ^ nl 

■^0,nZ r\ 


i+\i-n 


2 pSii/ini 'J'2 _ 

^ ) ^l,nl — 2 


nl 


nl 


We will average over all spin-scalar phases ijjni of the transmission amplitudes as there 
are usually many scattering channels in an area comparable with the superconducting 


coherence length squared. This hlters out all the terms in Eqs. (64a)-(646) where these 
scalar scattering phases cancel. 

For a magnetic system, in linear order in Tni and dnn' we obtain 
2(7 W 

r" = "" = 5,E„, [(To.„ii + Ti,„,K)G(T;,„,i + t;,„,k),g] 

-6,E„<«[k'.G]. (68) 

where Qq = e^/h is the conductance quantum. After multiplying out we obtain the set 
of boundary conditions 

2/W = [g^G + g^^{TG}+g^kGk-ig^k',G]^ (69a) 

with 


6 “ =e,E„,r„, 1 + ,/i-p, 


nl 


(696) 


= ^.EnZ^nZ (1 - VI - ^^Z j (69c) 

^MR ^ g^Y.nlTnlVnl, 6^“^ = 2^,En^nn (69d) 

For K = k' and the assumption of a channel-diagonal scattering matrix {n = 1) this also 
provides the derivation of the boundary conditions used for Ref. IST]. We now proceed 
to the second order terms: 




2/(2) = - 2/4 + gt [k'Gk\ G] „ + 1 
ME, =gl (G O Gk'+ k'G O G) +^5 Go [k',G] 

= gl {kGk o Git' + it'G O ztGzt) +glkGo [E, G] 


(70 a) 


x>x 


M 


M 


x,x 

MR 

x,x 




= g^^ {{TG} o Gk' + k'G O {k,G}) + {k,Go [E,G]} 
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where I 4 , denotes a cumbersome expression in fourth order of the transmission 
amplitudes, which we do not write down here explicitly (see 
used the abbreviations 


Appendix F). We have 




1 + 




~ Inl 


'dnn'Tnl I 1 






a 


MR 


= -^Gq^nl'^rinTnl'Pnh 




if 

./ U, 


2 

nn' 


(706) 

(70c) 

(70 (i) 

Note that 


and are defined as and with dnn replaced by 

ipnn' and do not appear in these expressions, in accordance with the intuitive notion 
that scalar scattering phases should drop out in the quasiclassical limit, which operates 
with envelope functions only. 

The case for only channel-conserving scattering (channel-diagonal problem) follows 
by taking in Eqs. {69b)-{69d) and (706)-(70d) only the terms with n = 1. All other 
formulas (69a), (70a) remain unchanged. This case is treated in Ref. [33] to linear order 


in Vnn, and our formulas reduce to these results for the considered limit. Note that for 
this case all spin-scalar phases cancel automatically and no averaging procedure over 
these phases is necessary. 


5. Application for diffusive superconductor/half metal heterostructure 

The problem of a superconductor in proximity contact with a half-metallic ferromagnet 
has been studied within the frameworks of Eilenberger equations [20l [11], [121 EO] [32] 
[67] l68] l69] , Bogoliubov-de Gennes equations [70] [TU [72] [73] , recursive Green function 
methods m, circuit theory na, within a magnon assisted tunneling model m, and in 
the quantum limit im. Various experiments on superconductor/half-metal devices have 
been reported, both for layered systems involving high-temperature superconductors 
[TSl [73] EOl ED and in diffusive structures involving conventional superconductors 

An important consequence of the new boundary conditions 


in Eq. (69a) is that half-metals can now be incorporated in the Usadel equation. 


appropriate to describe the second class of experiments mentioned above, whereas 
there previously existed no suitable boundary conditions to do so. Gonsider first a 
superconductor/half-metal bilayer with the interface located at x = 0 (see Fig. [^. 

The superconductor is assumed to have a thickness well exceeding the 
superconducting coherence length. Our expansion parameters are the spin-dependent 
reflection phase shifts at the superconducting side of the interface, dn,, and the tunneling 
probabilies Tni- For calculating triplet components in the half-metal it is sufficient to 
expand the solution for the Green function in the superconductor up to linear order, and 
the solution for the Green function in the half-metal up to quadratic order. The zeroth 
order term in the superconductor is pure spin-singlet, and the first order term pure spin- 
triplet. Thus, up to ini ending first order we can assume a bulk singlet order parameter. 
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Superconductor 



Half-metal 


Figure 2. A superconductor/half-metal bilayer with a magnetically inhomogeneous 
barrier region. The magnetization direction associated with the spin-dependent phase- 
shifts occurring on the superconducting side (described by the matrix k') does not 
in general align with the magnetization direction associated with the transmission of 
quasiparticles across the barrier (described by the matrix k). 


not affected by the interface scattering (corrections to the singlet order parameter arise 
only in second order in dn, and Tni)- For fntnre reference, we dehne the qnantities 
c = cosh(z/) = —1^, s = sinh(i/) = with v = atanh(|A|/i?), and 

denote the SC phase as 6 . We hnd for the triplet component in the superconductor 

, G'^cs 


Ftoi 


X) = i 


ascM 




ana. 


y 


(71) 


with the normal-state conductivity cxsc = 2e^iVsc-Dgc in the superconductor (Age and 
Dsc are the normal-state density of states per spin projection at the Fermi level and 
the diffusion constant, respectively), contact area A, and q = ^^/2D/hDsc- 

In the half-metal (width d), only spin-j' particles have a non-zero density of 
states at the Fermi level. In the spirit of quasiclassical theory of superconductivity, 
a strong exchange held is incorporated not in the transport equation, but directly in the 
band structure which is integrated out at the quasiclassical level [niE9], leaving only 
parameters like diffusion constant, and normal state density of states at the Fermi level 
for each itinerant spin band. For transport in a half-metallic ferromagnet, this means 
one must just include one spin-band with diffusion constant Fhm in the Usadel equation. 
Thus, only the elements and F|-|- exist in the Green function G of the half-metal. As 
we expand in the tunneling probability, we can (for energies well exceeding the Thouless 
energy h^nivi/d^ of the half-metal) use the linearized Usadel equation. 


2iFF|-'|- — 0. 


(72) 


Since there is only one anomalous Green function in the half-metal, we omit the 
spin indices for brevity of notation and dehne F = F^-|-. The general solution is 
F{x) = with A, B being complex coefficients to be determined from 

the boundary conditions, and k = 2iE/hD}s~A- At the vacuum edge of the half-metal 
{x = d), we have d^F = 0. At the interface between the superconductor and half-metal, 
the boundary conditions for F from the half-metallic side is obtained from Eqs. (69a)- 
{70d) with Vni = 1- Note that for Vni = 1, we have G^ = G^ = G^^ = Gx well as 
G'^ = G^ = We hnd that in order to obtain a non-vanishing proximity ehect, it 

is necessary that the magnetization direction associated with transmission across the 
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barrier (k) and spin-dependent phase-shifts picked up on the superconducting side of the 
interface (k') are different. We set k = dz since the barrier magnetization determining 
the transmission properties is expected to be dominated by the half-metal magnetization 
which points in the ; 2 -direction. The boundary condition for F aX x = 0 reads: 

anMAdz;F = 2iQ^cse'^{m^-im), Gv = Wx + a ( 73 ) 

(Jsc-TQ' 

with the normal-state conductivity cthm = e^iVHM-DnM in the half-metal is the 

normal-state density of states at the Fermi level), and the conductance Gv contains two 
terms: 2 Gx which is proportional to and a second term containing Qf^G^ 

which is proportional to (X)zl^zz)(X)nZ'7^^')- Moreover, rnf^ and mfy are the normalized 
components of a possible misaligned barrier moment compared to the magnetization of 
the half-metal. We have taken this into account by writing: 


K = 


+ mfy 


-e'^im' 


Inserting the general solution of F into the boundary conditions, one arrives at the hnal 
result for the proximity-induced superconducting correlations F in the half-metal: 

2 cosh[ik{x - d)] Gvcs , 

Fix) = - . ■ , , „ - [m., - imD. (75) 

^ ’ smh(iA;d) anmAk ^ ' 

This is the first time the Usadel equation has been used to describe the proximity effect 

in a superconductor/half-metal structure. Several observations can be made from the 

above expression. For small E the energy factors c oc. E in the numerator and k"^ (x E 

in the denominator cancel, such that the proximity-effect, if present, happens even at 

E = 0. The proximity-effect is seen to be non-zero only if spin-dependent scattering 

phases at the superconducting side of the interface are present, and at the same time 

their quantization axis k' is misaligned with that of the transmission amplitudes, k. 

The reason for this is that phase-shifts on the half-metallic side are irrelevant on the 

quasielassical level, because they are spin-scalar (only spin-)' particles have a hnite 

density of states there). On the other hand, the phase-shifts on the superconducting 

side have two consequences: they are responsible for an S'-m = 0 spin-triplet component 

on that side of the interface (where S is the spin vector of the Cooper pair), and 

they affect also transmission amplitudes. As a consequence, during transmission the 

quantization axis k' can be rotated into the 5*^ = ±1 spin triplet components which 

are allowed to exist in the half-metal if spin-flip processes exist at the interface {e.g. 

due to some misaligned interface moments). This is exactly the reason for why F also 

depends on and m/y whereas it is independent on the barrier moment m//. only a 

barrier moment with a component perpendicular to the magnetization of the half-metal 

can create spin-flip processes which rotate the S ■ m = 0 into the = ±1 components, 

and thus F also vanishes if m' = mi, = 0. 

Another important observation that can be made from the above expression is that 
a misaligned barrier moment effectively renormalizes the superconducting phase. Using 
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spherical coordinates, we may write rnf^ — imfy = sin 0 'e“^!^' where (p' is the azimuthal 
angle describing the orientation of the barrier moment in the xy-plane. Thus, the 
effective phase becomes 6 ^ 6 — ip'. To see what consequence this has in terms of 
measurable quantities, we proceed to consider a Josephson junction with a half-metal 
by replacing the vacuum boundary condition ai x = d with another superconductor. 
Solving for the anomalous Green function F in the same way as above, we may compute 


the supercurrent flowing through the system via the formula [see Eq. (15)]: 
cNumDhm-A. 


I = 


dETr{r3(GHM53;GHM) }• 


(76) 


Here, Tr denotes a trace over 2x2 Nambu-Gor’kov space. After some calculations, one 
arrives at the result: 


/ = Jo sin 0'^ sin sin( 6 'R — 9l + ip' 


J-L 


A) 


(77) 


where Jq is a lengthy expression depending on parameters such as the width d of the 
half-metal and the temperature T (and which vanishes unless and are non-zero). 
To be general, we have allowed the spin-dependent phase-shifts for each superconductor 
and the barrier moment at each interface to be different, indicated by the notation ’L’ 
and ’i?’ for left and right. We hnd that Jq is negative, giving rise to a vr-Josephson 
junction behavior for the case of = ip'^. 
ballistic case result of Refs. 


Expression (77) is consistent with the 


and shows how a hnite supercurrent will 
appear in a ring geometry even in the absence of any superconducting phase difference, 
~ = 0 , if the barrier moments are misaligned in the plane perpendicular to the 

junction, ip'^ — ip'^ 7 ^ 0 . A similar effect was also reported via circuit theory for a 
diffusive system eg. however not due to spin-dependent scattering phase shifts but due 
to some “leakage terms”. Within our formalism, we thus obtain a so-called 0o Josephson 
junction behavior [901 ED ED |93l [9l] with 0o = (vr -|- ip'^ — (pymod(27r). 

The above framework can be readily generalized to cover strongly spin-polarized 
ferromagnets building on the same idea as Ref. im. For a sufficiently large spin¬ 
splitting, the f- and j,-conduction bands can be treated separately in the bulk with a 
separate Usadel equation for and F^. These would then only couple via interface 
scattering and the strong exchange held would only enter by having different normal- 
state density of states A^.^, and diffusion coefficients D^, of the spin-bands in each 
separate Usadel equation. 


6. Conclusions 


We have derived new sets of boundary conditions for Usadel theory of superconductivity, 
appropriate for spin-polarized interfaces. We present a general solution of the problem 
appropriate for arbitrary transmission, spin-polarization, and spin-dependent scattering 
phases. The explicit equations for the most general set of boundary conditions are 
given in Eqs. ([M|)-(|345), (IS^-diO^, and (|4^-(|4^. With the solution of this long¬ 


standing problem we anticipate a multitude of practical implementations in future to 
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tackle superconducting systems that involve strongly spin-polarized materials. We have 
applied the general set of equations to various special cases important for practical 
use. We derived boundary conditions for an interface between a superconductor and a 
ferromagnetic insulator valid for arbitrary spin dependent scattering phases, Eq. (49). 
This extends previous work of Ref. [SS], which was restricted to small scattering phases. 


Using an exact series expansion of the general set of boundary conditions, Eqs. (50)-(55), 
we have obtained a perturbation series for the boundary conditions appropriate for such 
an interface, which allows for channel off-diagonal scattering and channel-dependent 
spin quantization axes, Eqs. {57a)-(57b). For the tunneling limit, we have presented a 
new set of boundary conditions appropriate for arbitrary spin polarization, non-trivial 
spin texture across the interface, and allowing for channel off-diagonal scattering, Eqs. 
(64a)-(64&). Neither of these three ranges of validity has been covered previously. As an 
application we then proceed to give a theoretical foundation of the boundary conditions 
used in Refs. [571 EH [59], Eqs. ( |69a[ )-( [M^ , which we have generalized for channel 
off-diagonal scattering and non-trivial spin texture across the interface. One central 
result of the application of our formalism is the extension of these relations to second 
order, including the important mixing terms between transmission and spin-dependent 
scattering phases. These terms, Eqs. (70a)-(70d) generalize the corresponding terms 
from Ref. [5^ to arbitrary spin polarization, possible nontrivial spin-texture across the 


interface, and channel off-diagonal scattering. We have demonstrated the application 
of the new set of boundary conditions by treating a diffusive superconductor/half- 
metal proximity junction and a diffusive superconductor/half-metal/superconductor 
Josephson junction. In the latter case we found a realization of a 0o-junction. We 
are confident that our boundary conditions will advance the field of superconducting 
spintronics considerably. 
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We perform a singular value decomposition of the reflection and transmission matrices 
(with dimensions n x n for ^n, m x m for S 22 , n x m for ^ 12 , and m x n for ^ 21 ) 


S 


^11 ^12 

*§21 —S22 


URV^ WTZ^ 

wfz^ -URV^ 


(A.l) 


\ / \ / 

Here U, V, W, Z, U, V, W, Z are unitary matrices, and the R, T, R, T contain the real and 
non-negative singular values in the main diagonal and are zero otherwise, he., T'l' = 
and T'l' = T^, R^ = R and R^ = R. We assume that the singular values are sorted from 
smallest to largest in R and R, and from largest to smallest in T and T. We introduce 
the unitary matrices , 'h = ZW <h = W^^IJ, and ^ = Z^'V. In terms of those, 

unitarity of the matrix S requires that (we denote for simplicity the unit matrices Inxn 
and Imxm with the same symbol 1; the dimension is clear from the context) 


(1 - R^) = ('A. 2 ) 

(1 - R^) = = ¥t^T^. (A.3) 


We see that 1 — R^ and 1 — R^ contain the eigenvalues of the Hermitian matrices on 
the right hand sides of the equations, which requires that these eigenvalues coincide 
with the values in the diagonal matrices TT\ T'I'T, TT\ and T'^T, respectively. Thus, 
with the sorting arrangement done above, the relations (1 — i?^) = TT'^ = T'^T and 
(1 — FR) = TT'I' = T^T hold. Because all singular values of T are real, this means that 
f = T\ R= Vl - TTt = Vl - f tf, R = Vl - TtT = Vl - ff t, and RT^ = TR, 
RT = T^R. Furthermore, the unitary matrices $ and T commute with R and the 
unitary matrices <F and T commute with R. In particular, those matrices are block 
diagonal, with block sizes given by the degeneracy of the singular values in R and R, 
respectively. The remaining unitarity requirements, using the above flndings, reduce to 


$Tt(T.R) = (T.R)T<ht (A.4) 

T$t(/2T) = (A.5) 


That means that for the blocks corresponding to non-zero reflection singular values 
the above two equations lead to the one condition <|>tTT = T'l’T'h. If there are no 
zero-reflection singular values then, remembering that <h commutes with R and T with 
R, 


S = 


0 

0 f/Tt 


0 

0 TVt 


(A.6) 


X 


R T 
Tt -R 

' X \ z X 

The blocks with zero-reflection singular values can be treated separately, and it is easily 
seen that the singular value decomposition of the scattering matrix has the general form 

Vl - TTt T 
Tt -Vl - TtT 

’ X ^ z X \ z 

with unitary matrices U, U, V, and V. The decomposition is not unique. 


S = 




(A.7) 
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An important featnre of the above representation is that the center matrix is Hermitian. 
If we only reqnire this property of the central part, but not necessarily diagonality of 
the m X n matrix T, then we can hnd an entire class of transformations that keep this 
property. We dehne TZDTZ'^ = T with unitary matrices TZ and TZ. Then 


S = 


0 W Vl - DDt D \ 0 A 

0 UiZ ) \ Dt -yJi-DW ] I 0 'RfV^ ] 

(B.l) 


where D is now an n x m matrix that is not necessarily diagonal anymore. Consider 
now some special cases. First, we chose TZ = V\ TZ = VC Then 

UV^ ° ^ Vl - C'C'^ C' 

0 WVt J i c'^ 


S = 


(B.2) 


C'^C 

’ ^ \ ^ 
with C = VTVt gives a polar decomposition of the reflection parts of the scattering 
matrix S. Similarly, TZ = U.^^ Iz = leads to 

Vi - CCt c 

-vT 


s = 


Ct 


ac 


UV 0 
0 WVt 


(B.3) 


X 


X 


with C = UTU^ = WV^C"(WVl)l. We can also chose a decomposition in the form 

Vl - C"C"^ C" 


s = 






(B.4) 


with C" = VTU\ or other decompositions. 

These decompositions are unique when there are no zero-reflection singular values. 
This means, that under the conditions of no zero-reflection channels and are 
uniquely defined, as the matrices C and D are. The unique unitary matrices UV"' and 


Wyl are the surface scattering matrices S and S, Eq. (19). 


Appendix C. Parameterization of scattering matrix 


We now turn to a useful parameterization of the transmission matrix C. We note that 
with the definition 


we obtain 


C = (1 + ^2t 


Vi - ccf c 
cl -Vi-cfc 


X 



(C.l) 


(C.2) 
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with 


f = (l + tt^) ^ (l — tt'^) 
r = (l + th) ^ (l — th) 
d = (l + tt^) ^ 2t. 


(C.3) 

(C.4) 

(C.5) 


To connect with the main text, see equations (20)-(21). Furthermore, if t = uOv"^ is a 
singular decomposition for t, then C = u[{l+9‘^)~^29\v'^ is a singular decomposition of C. 
Conversely, if C = u5v^ is a singular decomposition for C, then t = m[( 1 — -y/l — 6‘^)/5]v^ 
is a singular decomposition for t. If 0 < 0 < 1 then 0 < 5 < 1 and vice versa. Thus, the 
parameterization in terms of t is equivalent to that in terms of C. 


Appendix D. Expansion to third order of expression (56) 


To third order we obtain from Eq. (56) 
= 1 [k, G] 


and 


X(2) = _ _^K^G]+i[k,{G + ivri) o o G] 

X(3) = _^_^KSG]-^[kUG + iTTl) o (G')(') o G] 

+ i[A,(G + i7ri)o(G')^^)oG] 

+ 1 [k, {G + ivri) o oGo (G')^i) o G] 


(«')“’ = -k X.G] . (G')'"' = - A [K ? G] . 


27r^ dvr^ 

This can be simplified further noting 

60 ( 6 ')^ =_(g')(i)oG, 

O (G')^i) = - {(G")(2)^(5| ^ 

(G + ivri) o G = i7r(G + i7ri), 

271^ [k ^ (G')^^^] = -i[k G] , 

4^2 [A-, (G')'">] = - [A ? G] , 


yielding Eq. (57a) of the main text. 


(D.l) 

(D.2) 


(D.3) 

(D.4) 


(D.5) 

(D.6) 

(D.7) 

(D,8) 

(D.9) 


Appendix E. Expansion of matrix current for finite transmission 


From section |4.3| we obtain the following expressions to second order in the spin 
dependent reflection phases and in the transmission probability: 

i(‘) = -27n[(“,jS]7 + i[A,G], 

= -2X1 [A, jS]f + > [A. (A)'*’] - ^ [A ? G] , 


(E.l) 

(E.2) 
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and 
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P.sSll'’ = [foQflG ^ 

P.sSll'’ = [G..G]„ 

P.sSlf = +[Gi.(9S)“’]. 

L ^ Jo 

+ G10G0 [Gi,G]^oGoGi 

+ i ([[*-, Gi] , G]^ + [t„ [£, G] T„', g] J , 


(§•)<« = (G + ijri) o ([Gi, G] ^ + (G')<‘' o G) 
(a)" = (G+i!rl)o(G)“oG 
( 90 °)'“’ = (G - iiri) o (G")<‘> o G 


(E.3) 

(E.4) 


(E.5) 

(E.6) 

(E.7) 

(E,8) 


with from Eq. (D.4). Collecting everything together, we obtain the 

resnlt shown in Eqs. (63a)-(636) of the main text. 


Appendix F. Term of second order in transmission probability 

For completeness we present here the expression of order Tfi- 
h = glGoGoG + gl kGk o G o kGk 

+ gf^ (G o G o kGk + kGk o G o G) + gf^' {k,G} o G o {k,G} 

+ gT" {GoGo{k,G} + {k,G}oGoG) 

+ gf^' {kGkoG o {k,G} + {;t,G} o G o kGk) (F.l) 

with 

g^ = Pnln'v'TnlTn'l' ^ + \/^ ~ '^n'l' ) (^-2) 

nln'V 

si = gS, I] (i - y'l - VI ) (i - fi-vi ,,,) (f.s) 

nln'V 

E Pnln-l-T„,T„-V (l + f-K ) (i - di - ril') (E-'i) 

nln'V 

gV^ = -^gq Pnlnn'TnlTn'v'PnfPn'V (F-5) 

nln'V 

Pnln'VrnlTn'vVn'V (l + \f^^l ) (F.6) 

nln'V 

Y 1 Pnln'vTnlTn'l'Vnn' (l " ) (F.7) 

nln'V 

with Pnin'i' = dnn' + dw — dnu'div, arising from averaging over the typical phase 
factor of spin-scalar transmission phases. The channel-diagonal case 

follows from setting n = I = n' = I' and Pnnnn = 1- 
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